KOLAKOSKI-(3,l) IS A (DEFORMED) MODEL SET 



MICHAEL BAAKE AND BERND SING 

Abstract. Unlike the (classical) Kolakoski sequence on the alphabet {1, 2}, its analogue on 
{1,3} can be related to a primitive substitution rule. Using this connection, we prove that 
the corresponding bi-infinite fixed point is a regular generic model set and thus has a pure 
point diffraction spectrum. The Kolakoski- (3, 1) sequence is then obtained as a deformation, 
without losing the pure point diffraction property. 



1. Introduction 

A one-sided infinite sequence u over the alphabet A = {1, 2} is called a (classical) Kolakoski 
sequence (named after W. Kolakoski who introduced it in 1965, see |H_), if it equals the 
sequence defined by its run lengths, e.g.: 

u = ^22, > U y 2 ^_L^ 22 ^1^ ^22, JA^ 2 ^11^ ... 

^ ' 2 2 1 1 2 1 2 2 1 2 ... = u. 

Here, a run is a maximal subword consisting of identical letters. The sequence u' = loo is the 
only other sequence which has this property. 

One way to obtain lu of (JJ) is by starting with 2 as a seed and iterating the two substitutions 



1 i — ► 2 1 i — ► 1 

2 ^ 22 and ai : 2 h+ 11, 



alternatingly, i.e., <r substitutes letters on even positions and letters on odd positions (we 
begin counting at 0): 

2 i — ^ 22 i — ^ 2211 i-> 221121 ■-> 221121221 h+ . . . 

Clearly, the iterates converge to the Kolakoski sequence oj (in the obvious product topology), 
and ijj is the unique (one-sided) fixed point of this iteration. 

One can generalize this by choosing a different alphabet A = {p, q} (we are only looking 
at alphabets with card(.A) = 2), e.g., .A = {1, 3}, which is the main focus of this paper. Such 
a (generalized) Kolakoski sequence, which is also equal to the sequence of its run lengths, can 
be obtained by iterating the two substitutions 

q i— ► p q . q i— > q q 

<7 n : „ and Ci : „ 

p I — > pr x p I — > q F 

alternatingly. Here, the starting letter of the sequence is p. We will call such a sequence 
Kolakoski- (p, q) sequence, or Kol(p, q) for short. The classical Kolakoski sequence to of Q is 
therefore denoted by Kol(2, 1) (and to' by Kol(l, 2)). 

While little is known about the classical Kolakoski sequence (see |15j). and the same holds 
for all Kol(p, q) with p odd and q even or vice versa (see ^ ne situation is more favourable 
if p and q are either both even or both odd. If both are even, one can rewrite the substitution 
as a substitution of constant length by building blocks of 4 letters (see [^E21)- Spectral 
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properties can then be deduced by a criterion of Dekking [T^]. The case where both symbols 
are odd will be studied in this paper exemplarily on Kol(3, 1). 

It is our aim to determine structure and order of the sequence Kol(3, 1). This will require 
two steps: First, we relate it to a unimodular substitution of Pisot type and prove that 
the corresponding aperiodic point set is a regular generic model set. Second, we relate this 
back to the original Kol(3, 1) by a deformation. Here, the first step is a concrete example of 
the general conjecture that all unimodular substitutions of Pisot type are regular model sets 
(however, not always generic). This general conjecture cannot be proved by an immediate 
application of our strategy, but we hope that our method sheds new light on it. 

Remark: Every Kol(p, q) can uniquely be extended to a bi-infinite (or two-sided) sequence. 
The one-sided sequence (to the right) is Kol(p, q) as explained above. The added part to the 
left is a reversed copy of Kol(q,p), e.g., in the case of the classical Kolakoski sequence of (0), 
this reads as 

. . . 11221221211221|22112122122112 

where "|" denotes the seamline between the one-sided sequences. Note that, if q = 1 (or 
p = 1), the bi-infinite sequence is mirror symmetric around the first position to the left 
(right) of the seamline. The bi-infinite sequence equals the sequence of its run lengths, if 
counting is begun at the seamline. Alternatively, one can get such a bi-infinite sequence by 
starting with q\p and applying the two substitutions to get a 1 (q)\a (p) in the first step and 
so forth. This also implies that Kol(p, q) and Kol(q,p) will have the same spectral properties, 
and it suffices to study one of them. 

2. KOL(3, 1) AS SUBSTITUTION 

If both letters are odd numbers, one can build blocks of 2 letters and obtain an (ordinary) 
substitution. Setting 1 A = 33, B = 31 and C = 11 in the case of Kol(3, 1), this substitution 
a and its substitution matrix M (sometimes called incidence matrix of the substitution) are 
given by 

A i — ^ ABC / 1 1 1 \ 

(2) o~ : B i — ^ AB and M = 1 1 , 

C ■-> B \ 1 / 

where the entry Mjj is the number of occurrences of j in a(i) (i,j £ {A, B,C}; sometimes 
the transposed matrix is used). A bi-infinite fixed point can be obtained as follows: 

(3) B\A i — ► AB\ABC i-> ABCAB\ABCABB ■-> ... 
This corresponds to 

(4) ... 3331113331|333111333131 . . . 

which is the unique bi-infinite Kol(3, 1) according to our above convention. The matrix M is 
primitive because M 3 has positive entries only. The characteristic polynomial P(x) of M is 

(5) P(x) := det(x 1 — M) = x 3 — 2x 2 — 1, 

^That Kol(3, 1) can be related to a substitution is well-known, e.g., in |14| . a substitution over an alphabet 
with four letters is given, while I33| uses the same substitution with three letters as we do. We thank the 
referee for pointing this last reference out to us. 
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which is irreducible over Z (there is no solution mod 3) and over Q (every rational algebraic 
integer is an integer). The discriminant D of P(x) is D = so P(x) has one real root a 
and two complex conjugate roots (3 and 0. One gets 

2.21 « a > 1 > \/3\ » 0.67 > 0, 

wherefore a is a Pisot-Vijayaraghavan number (i.e., an algebraic integer greater than 1 whose 
algebraic conjugates are all less than 1 in modulus), and a is a substitution of Pisot type. Since 
det(iW) = 1, the roots a, (3 and f3 are also algebraic units, and the associated substitution 
is said to be unimodular. Note that Re(/3) = 1 — %. If necessary, we will choose (5 such 
that Im(/3) > in the following calculations (the other possibility only leads to overall minus 
signs) . 

There is a natural geometric representation of such a substitution by inflation, compare |24j . 
Here, one associates bond lengths (or intervals) £ A , £ B and £ c to each letter. These bond 
lengths are given by the components of the right eigenvector which belongs to the (real) 
eigenvalue a and is unique (up to normalization) by the Perron-Frobenius theorem. The 
normalization can be chosen so that 

£ A = a 2 - a « 2.66, £ B = a w 2.21 and £ c = 1. 

Inflating the bond lengths by a factor of a and dividing them into original intervals just 
corresponds to the substitution (because a ■ £ A = £ A + £ B + £ c , etc.). We will denote this 
realization of the bi-infinite fixed point with natural bond lengths (respectively the point set 
associated with this realization where we mark the left endpoints of the intervals by their 
name) by £Kol(3, 1), reserving "Kol(3, 1)" for the case of unit (or integer) bond lengths. 

On the other hand, the frequencies p A , p B and p c of the letters in the infinite sequence 
are given by the components of the left eigenvector of M to the eigenvalue a. This gives 

p A = ^(-a 2 + 3a-l) «0.38, p B = a 2 - 2 a ss 0.45, p c = - (-a 2 + a + 3) « 0.17, 

with p A + p B + p c = 1. Therefore, the average bond length £ in the geometric representation 
is 

(6) £ = p A .£ A + p B .£ B + p c .£ c = l - (-a 2 + a + 7) « 2.17, 

and the frequencies of 3s and Is in Kol(3, 1) can easily be calculated to be p 3 = — 1) ~ 0.60 
and p x = |(-a + 3) » 0.40. 

Remark: In the case where p and q are odd (positive) integers, one gets unimodular sub- 
stitutions of Pisot type iff p = q ± 2. More generally, one gets substitutions of Pisot type iff 
2 ■ (p + q) > (p — q) 2 holds. Otherwise, all the eigenvalues are greater than 1 in modulus, 
see [3*T] . 
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3. Model Set and IFS 

A model set A(Q) (or cut-and-project set) in physical space M. d is denned within the following 
general cut-and-project scheme [2H E] 



TT 



x H — ► H 



(7) 



\ U / 

1-1 dense 

r 



where the internal space H is a locally compact Abelian group, and r C M. d x iJ is a lattice, 
i.e., a co-compact discrete subgroup of M d x il. The projection w int (r) is assumed to be dense 
in internal space, and the projection n into physical space has to be one-to-one on r. The 
model set A(Q) is 

A(Sl) = {tt(x) \ xeT, vr int (x) G n} C M d , 

where the window C H is a relatively compact set with non-empty interior. If we set 
L = Tt(r) C M d , we can define, for x £ L, the star map * : L ^ by x* = 7r int o (Tr^)^ 1 (x), 
see [SJ- So we have r = {(x,x*) \ x E L} and L* = 7r int (.T). If the boundary 90 of the 
window has vanishing Haar measure in H, we call yl(fi) a regular model set. If, in addition, 
7r int(-^ n ) H dQ = 0, the model set is called non-singular or generic. 

Every model set is also a Delone set (or Delaunay set), i.e., it is both uniformly discrete 2 
and relatively dense 3 . A Delone set X is a Meyer set, if also X — X is a Delone set. Every 
model set is a Meyer set, see [25] . 

We will now construct a model set A(JY) and - in a first step - show that this model set 
differs from SKol(3, 1) at most on positions of density 0. By Galois conjugation (see pHITT] b 
which here corresponds to the star map as we will see, we find a lattice 

(8) r = Z • v A + Z • v B + Z • v c C R x C ~ M 3 

where 

Q: 

u A = I Re(/^ - 0) | , « B = | Re(/3) | and v c 





Im(/3) 

The projection ir (i.e., the projection on the first coordinate) is injective on r because Q(a) 
is a Q-vector space of dimension 3 with (Q-) linearly independent elements 1, a and a 2 . Also, 
ir(r) = Z[a] is dense. To see that 7r int (.r) is dense, we note that ir int (r) = Z[/3] and that 1 
and (3 are linearly independent. So, (3 n and /3 n+1 are also linearly independent for all n G N, 
and their Z-span forms a two-dimensional lattice in C, which is a uniformly discrete subset 
of Z[/3]. Since \f3\ < 1, one can choose, for every e > 0, an n, such that there is a lattice point 
(of the lattice Z • /3 n + Z • /3 n+1 ) in every ball of radius e, so Z[/3] is dense in C. Note, that 
7Ti nt is also injective on r (this can be seen from Re(/3) = 1 — f and Re(/3 2 ) = 2 — ^-). So we 
have established: 



2 A set A is uniformly discrete if 3r > s.t. every open ball of radius r contains at most one point of A. 
3 A set A is relatively dense if 3R > s.t. every closed ball of radius R contains at least one point of A. 
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Proposition 1. With C ~ M?, r of @ and the natural projections it and ir int , we obtain the 
following cut-and-project scheme: 



(9) 



7T 



dense ^ ' dense 

1-1 1-1 



' dc 

r 



Furthermore, we have 

7r(r)=Z[a] and vr int (r) = Z[/3], 
where a is the real root of © and (3 one of the complex conjugate ones. □ 

In order to describe £Kol(3, 1), the main task is now to determine the appropriate windows 
Q, A , Q B and £l c (one for each letter; f2 = Q A L)Q B U£l c ). For these windows, the substitution 
rule a of © induces the following iterated function system (IFS for short) in internal space , 
cf. EH: 

n A = pn A u pn B 

(10) n B = (3n A + (3 2 -(3 u pQ B + f3 2 -p U (3Vt c 

n c = pn A + p 2 . 

This IFS is obtained as follows: We denote by A A the subset of £Kol(3, 1) of left endpoints 
of intervals of type A (of length £ A ), and similar for A B and A c (we have XKol(3, 1) = 
A A \jA B UA c , where U denotes disjoint union). Then, the substitution a of (J2J) induces the 
following equations for these Delone sets in M: 



c 



A A = a A A U a A B 

(11) A B = aA A + a 2 -a U aA B + a 2 -a U a A, 

A c = aA A + a 2 . 

Applying the star map to these equations yields (|1()|) . In this sense, the iteration of the 
IFS (jl()j) in internal space corresponds to the iteration (|lljl in physical space (note that by 
Proposition n the star map is bijective on Z[a]). 

Setting Q AB = Q A U $7 B in (|1()|) , the system decouples and we remain with the simpler IFS 



(12) n AB = A(^) u f 2 (n AB ) u f 3 (n AB ) 

where 

(13) h(z)=0z, f 2 (z) = p 3 z + p 3 and f 3 (z) = (3 z + /3 2 - 0. 

The mappings f i : C \—* C are contractions (|/3| < 1), so that Hutchinson's theorem [2 
Section 3.1(3)] guarantees a unique compact solution of ()12|) . called the attractor of the IFS 
The sets fl A , tt B and £l c can be calculated from Sl AB as 

(14) o A = h(si AB ), n B = f 2 (n AB )uf 3 (n AB ), and n c = f A (n AB ), 

4 For later reference, we write: 

GZ3) n B = / 3 (n A ) u / 3 (n B j u A(n c ) 

«o = MOO, 
where f 1 and / 3 are defined as in 11131 . and f (z) — (5 z + f3 2 . 
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where / 4 (z) = (3 2 z + j3 2 . They are also compact sets in the plane. For the components of 
Q AB , see Figure E the windows 0, A , etc., are shown in Figure El Note that the decoupling of 
the IFS (i.e., the step from H10|) to l|12|)) lies at the heart of our argument and seems to be the 
reason that we cannot immediately generalize our method to other unimodular substitutions 
of Pisot type 5 , because no such decoupling emerges in general. 

The similarity dimension s of a set given by an IFS is the unique non-negative number s 
such that the contraction constants to the power of s add up to 1 (see JB])- For £l AB , this 
means 

|/?| 3s + 2|/3| s = 1 

with solution s = s(Cl AB ) = 2 (because a ■ \(3\ 2 = 1, the substitution is unimodular). The 
similarity dimension s of a set is connected to its Hausdorff dimension h by h < s where 
equality holds if the open set condition (OSC for short) is satisfied 16 . An IFS with mappings 
/j satisfies the OSC iff there exists a nonempty open set U such that f^U) n fj{U) = for 
i ^ j and f^U) C U for all i. It is easy to see that the corresponding self-similar set fi AB 
must be contained in the closure U so that the pieces fi(Q AB ) C f^(U) can intersect at their 
boundaries but cannot have interior points in common jSJ. If their boundaries do intersect, 
the IFS is called just touching. 

Proposition 2. The IFS of (| jjgj) for Q AB is just touching. 

Proof. To determine the boundary of Q AB , we choose special points P i in 0, AB , see Figure ^ 
(for illustration) and Tabled (for details) 6 . We first show how these points are determined. 
Demanding 

P2 = h(Pi), n = /iCP 2 ), Pa = /i(is)> and P 1 = f 3 (P 4 ) 
one gets the following fixed point equation 

(15) P x = f 3 o f\ o f x o / 3 (P X ) = /3 4 P 1 + 6/3 2 + 2 

for P l5 and similar results hold for P 2 , P 3 and P 4 . The unique solution of (|15[) is 

Px = ^(/3 2 - 3/3-1). 

Choosing 

^5 = ^1+^), ^Wsfo), P7 = fl(Pd, 

P 8 = f2(P 3 ), P 9 = fl(Pl), PlO = U P 2 + P 3 ) 

and setting r to be the inversion in the center P 5 (r : z \— * —z — f3) and k the one in the 
center P 10 (k : z i— ► — z — /3 + 1), one can verify the following equations: 

Pi = r(P 3 ) 

P 2 = r(P 4 ) = «(P 3 ) 

n = ^2 + PJ = 1(^6 + ^9) = 5(^7 + ^) 
P 6 = r(P 9 ) 

P7 = /2(A) = /3OP9) = r(P 8 ) 

P 8 = fi(P 6 ) = h( p 2) = r(P 8 ) = «(P 9 ) 

= 2(^8 + ^9)- 

^The substitution @ can be analyzed by the balanced pair algorithm as described in 1341 . This algorithm 
also confirms that it has pure point spectrum, but one does not get the model set property. 

6 We use the two dimensional geometry of the internal space here explicitly, and, instead of going into 
cumbersome notations and explanations, show some figures to clarify and assist the proofs. 
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Figure 1: Components fii^l^g) of the set 
and the points P 1 , . . . , P 10 used in the proof of 
Proposition [21 is the union of the three 

shaded areas. 



Point 


Coordinate 


Pi 


\{P 2 - 3/3-1) 


P2 


I (/? 2 - 3/3+1) 


P3 


i(-/3 2 + /3 + l) 


P, 


i(-/3 2 + /3-l) 


P 5 


-±/3 


Pe 


±(/3 2 -/?-l) 


Pi 


i(-/3 2 -/3-l) 


Ps 


±(/3 2 -/? + l) 


P9 


i(-/? 2 -/? + l) 


Pio 





Table 1: Coordinates of the points 
P 1 ,...,P w in Figure [TJ 



For the mappings, one finds 

/1 or = 7-0/3, / 2 °t = to/ 2 and / 3 ot = to/ 1 , 

showing that is inversion symmetric in the center P 5 , i.e., T(f2^g) = f2^ B . 

Denoting by [P 2 ,P 3 ] the "boundary" between P 2 and P 3 (the "right edge"), one finds 

[P2,P*] = /3 ^0/1(^3]) 

[P8,P Q ] = Aoro/joro/^Pj]) 
[^9^3] = Zioro/^Pj]) 
and therefore the following IFS for [P 2 , P 3 ] : 

(16) [P 2 ,P 3 ] = 5l ([P 2) P 3 ]) U 52 ([P 2 ,P 3 ]) U </ 3 ([P 2 ,P 3 ]), 
where 

(17) 9l (z) = -/3 2 z-/3, <? 2 (z) = (2/3 2 + l)z + /3 2 + l, and g 3 (z) = -/3 2 z - /3 2 . 

Of course, we have not shown yet that [P 2 ,P 3 ] really is (a piece of) the boundary of &a B , so 
we just define [P 2 ,P3] to be the unique compact solution of the IFS (|16|) . which is inversion 
symmetric in the center P 10 because 

J 1 ok = koj 3 , g 2 o k = K o g 2 and 5 3 ok = koj 1 . 

Also, we know that [P 2 , P 3 ] is connected since we can start the iteration with the straight line 
from P 2 to P 3 . In each iteration, the image remains a (piecewise smooth) path from P 2 to P 3 . 
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With the mappings f 1 and r, we get a boundary 

[P 2 ,P 3 ] U A([P 2 ,P 3 ]) U t([P 2 ,P 3 ] U to/^P^]) = 
[P 2 ,P 3 ] U [P 3 ,P 4 ] U [P^] U [P 15 P 2 ] 

around a simply connected open set U (we will prove in the next proposition that this bound- 
ary is non-self- intersecting) . Now one can show that only the boundaries of f i {U) intersect. 
Consider, for example, the region between f 2 (U) and f 3 (U). Then the boundary [P 7 ,P 8 ] on 

f2( u ) is g iven b y 

[P7,p*] = ^Wia^sD u / 2 ([p 2 ,p 3 ]), 

while the one on f 3 (U) is given by (taking orientation into account) 

[P 7 ,P 8 ] = f 3 o K og 2 ([P 2 ,P 3 ]) U f 3 og 3 ([P 2 ,P 3 ]). 
It is easy to verify that 

f 2 o T o f 1 = f 3 o ko g 2 and f 2 = f 3 o k o fif 3 . 

So the boundaries coincide. Similarly, one can check the region between fi(U) and f 2 (U) 
(note that P 7 and P 8 belong to all three sets f^U)), and that the boundary of U coincides 
with pieces of the boundaries of the f^U) - so the situation is as expected from Figure ^ 
Therefore, the IFS is just touching. Also, we now know that [P 2 ,P 3 ] is really a piece of the 
boundary of &ab- O 

Proposition 3. Let ^l^ B be the unique compact solution of the IFS (jj^jl . Then its boundary 
is non-self-intersecting. 

Proof. We specify a (closed) rhombus R (which surrounds the boundary [P 2 ,P 3 ] and the 
straight line from P 2 to P 3 ) such that its iteration in (|12j) will not leave R, i.e., such that 

(18) ( 9l (R)Ug 2 (R)Ug 3 (R))n(R 2 \R) = 0. 
Furthermore, we also require that 

(19) 9l (R) n g 3 (R) = 
and 

9l (R) n g 2 ( gi (R)) = forf/1 
r , m 92(R)ng 1 (g i (R)) = forf/3 

1 ' g 2 (R)ng 3 ( gi (R))=0 for i + 1 

g 3 (R)ng 2 ( 9i (R)) = for 3. 

Such a rhombus P exists, see Figure |^1 for a picture of such a rhombus that satisfies the 
conditions of (fTH)) and (fT§)) and Table [21 for the coordinates of its corners (of course, we could 
also use a shape different from a rhombus). This rhombus also satisfies ()2U(I . see Figure |3 

Here, 1)19(1 tells us that [P 2 ,Pg] and [Pg,P 3 ] do not have a point in common; similar state- 
ments apply for ()20() . Each iterate of the rhombus is associated to a corresponding iterate of 
the boundary [P 2 ,P 3 ] or the straight line from P 2 to P 3 (denoted by (P 2 ,P 3 )). We call two 
rhombi at the same iteration level neighbouring if their corresponding iteration of (P 2 ,P 3 ) 
have a common endpoint. We see in Figure that only neighbouring rhombi intersect at 
the second iteration level. We show that for any iteration level only neighbouring rhombi 
intersect. 
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Figure 2: A rhombus R 
with corners E 1 , . . . , P 4 , such 
that g^R) C R for all % 
(see (Pj|). The g^R) are 
shaded in gray. For the co- 
ordinates of E 1 ,...,E i , see 
Table E 




Figure 3: In addition to Fig- 
ure El the second iteration of 
the rhombus R is also shown, 
shaded in the dark gray. (|20|) 
is satisfied. 




Figure 4: The situation 
around the "joint" P 8 . We 
have an inversion symmetry 
in the center P 14 . For the co- 
ordinates, see Table |2j 



We have verified the assertion for the first and second iteration level and proceed induc- 
tively. Since g 1: g 2 and 53 are affine, we get the third iteration level as follows: The associated 
rhombi between P 2 and P 8 are a scaled down (by g x ) version of those of the second level, 
therefore the assertion holds for them. Similarly for the rhombi between P 8 and P 9 (by g 2 ) 
and between P 9 and P 3 (by g 3 ). So the only critical points remaining are the "joints" at 
P 8 and Pg. We show that at these points also only neighbouring rhombi intersect and for 
this, we make use of the self-similar structure of the boundary, see Figure |IJ The boundary 
[P 12 , Pu] is inversion symmetric in the center P 14 . This is clear for [P 13 , P 8 ] = ^(^([P^j -P3D) 
(and therefore P 14 = #1(53(^10)))- But [t also holds for [ p i2> p i3] = Si (02(^2) ^3])) and 
[P 8 ,P 11 ] = 02(51 ( K ([-^2> -^3])))- So, since the assertion holds around P 13 , it also holds around 
P 8 by symmetry. Similar arguments apply around P 9 . So the assertion holds for the third 
iteration level, i.e., for the third iteration level only neighbouring rhombi intersect. But the 
same argument applies to all further iteration levels. So the assertion is true, i.e., for a given 
iteration level only neighbouring rhombi intersect. Also note that each rhombus has two 
neighbouring rhombi (with the exception of the "starting" and "ending" rhombi at P 2 and 
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P 3 which only have one) and that there is no "rhombus loop", i.e., going from P 2 to P 3 we 
cross each rhombus only once. 

Now, suppose [i-^-Fy is self-intersecting. Then there exist points x,y 6 [P^-^y (x ^ y) 
such that they are connected in [P 2 , P 3 ] in two different ways, W 1 and W 2 (and we have a loop). 
We can choose points u £ W 1 and v G W 2 such that d(u,W 2 ) = mm{d(v,z) \ z £ W 2 } > 
(W 2 C [P 2 ,P 3 ] is compact) and d(v, W r 1 ) > 0. But then, u, v are in non-neighbouring rhombi 
for some iteration level N (and then for all iteration levels n > N), since the length of a 
rhombus of the iVth iteration level is at most \/3\ 2N • \E 2 — E 4 \. So, we get a "rhombus loop" 
for this iteration level by the rhombi which overlay W 1 and W 2 . This is a contradiction, 
therefore [^25-^3] is non-self- intersecting. 

From this single edge we proceed to all of the boundary. Here, critical are the "joints" 
P 2 , P 3 , etc., again, because we get the other three parts by an affine map of this edge 
(e-g-, [-P 3 ,-P 4 ] = /i([-P 2 ,P 3 ])) and opposite edges (i.e., [P 2 ,P 3 ] and [P 4 ,Pi]) do not overlap, 
cf. Figure El But at P 2 , an argument like the one at P 8 above applies, i.e., we have an 
inversion symmetry of part of the boundary in the center ^ (P 12 + -P 13 ) (and similar for the 
other "joints"). This extends our findings to the entire boundary. □ 

This also implies that the boundaries of Q A , £l B and £l c , respectively their union f2, are 
non-self-intersecting. Also, from the proof of the last proposition, we can deduce the following. 

Corollary 1. The point is an inner point of fi^^g) C £l A and —(3 is an inner point of 

f 3 (n AB )cn B . 

Proof. We again use the iteration of rhombi as in Proposition |3] to show that the two points 
are really inner points in the respective areas. For this, see Figure |SJ where the first iteration 
of the rhombi is used for all parts of the boundary. Clearly, the points 0, —(3 are inner 
points, which can easily be checked by a simple (though somewhat tedious) calculation of 
distances. □ 

Proposition 4. Let £l AB be the unique compact solution of the IFS (| H| . 

(i) £l AB is inversion symmetric in the center P 5 = —\f3. 

(ii) £l AB has Hausdorff dimension h(£l AB ) = 2. 

(iii) Q AB has positive (Hausdorff and Lebesgue) measure (area). 

(iv) The boundary dVL AB has vanishing (Lebesgue) measure. 

(v) There is a periodic tiling of the plane with £l AB as prototile. 



Proof, (i) See proof of Proposition |2 

(ii) Just touching implies the OSC, therefore h(Q AB ) = s(ft AB ) = 2. 

(iii) The OSC for £l AB (or any self-similar set with similarity dimension s) is equivalent to 
the positive Hausdorff measure condition fj, s (Q AB ) > 0, where /u s denotes the s-dimensional 
Hausdorff measure, see 8 and references therein. For Euclidean dimensions, Hausdorff and 
Lebesgue measure are connected by a nonzero multiplicative constant. 

The similarity dimension s = s(d£l AB ) of the boundary is the solution of (contraction 



iv 



constants given in (|17[)) 

2|/3| 2S + |/?| 3s ~ = l, 

which is s = — log(r)/ log(|/?|) ~ 1.22 (where r = |(1 + \/5) is the golden ratio; the previous 
equation is solved by \(3\~ s = r). Therefore, the statement follows from h(dQ AB ) < s(dQ AB ). 



(v) Because of the inversion symmetries r of $l AB and k of d£l AB from P 2 to P 3 (see proof of 
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Point 


Coordinate 


E x 


P 10 -z§(/3 2 -2/3) 


E 2 


P 2 + |(/5 2 -2/3) 




P 10 + z|(/3 2 -2/5) 


E, 


P 3 -|(/3 2 -2/3) 


P 8 


I(/? 2 -/? + l) 


Pu 


i(9/3 2 + /3 + 5) 


Pl2 


i(-3^ 2 - 3/3-1) 


Pl3 


i(5/3 2 -/3 + 3) 


Pu 


i(3/3 2 -/3 + 2) 



Table 2: Coordinates of the 
points used in Figures |21 
and|U 




Figure 5: The points and —j3 are represented 
as black dots. Also, the boundary is shown 
and the area of the rhombi (first iteration for 
all parts of the boundary) is shaded in gray. 
The boundary runs inside this shaded region, 
so and —(3 are inner points. 



Proposition EJ) , the "right edge" and the "left edge" differ only by a translation P 2 — Pi = 1, 
and, similarly, the "upper edge" and the "lower edge" differ by P 2 — P 3 = (3 2 — 2 (3. □ 

Proposition 5. Let £l A , £l B , £l c be the solution of the IFS <\l(Jty . and Q = £l A U Q B U 

Then Q is a compact set, homeomorphic to a disc, with positive area. The boundary dfl is 
a fractal of vanishing Lebesgue measure, which is non-self-intersecting. The set $7 admits a 
lattice tiling ofM?, where the lattice is spanned by P 2 — P 6 = —f3 + 1 and P 2 — P 3 = (3 2 — 2 0. 

Proof. It is clear from our construction that £1 is a compact set with simply connected interior. 
We have also seen that the boundary is connected and consists of finitely many pieces, each 
of which is obtained from a construction as used in the proof of Proposition |2I So, Q must 
be homeomorphic to a disc. The remaining statements follow directly from Propositions 
and|IJ because the mappings in (|14|) are affine and the just touching property also holds for 
f2 = £l A U f2 s U £l c . Since we also know the boundary of O (we have an IFS for every part 
of it), we can also verify the translation vectors by comparing the corresponding iterated 
function systems. Also, see Figure El for a depiction of these vectors. □ 



Corollary 2. A(£l) is a regular model set. 



□ 
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Ira 




Figure 6: Components of the set Q, the periodic tiling of the plane (C) with it and the 
corresponding translation vectors. 



We can calculate the volume 7 

(21) |P| = \deb(y A ,v B ,v c )\ = |Im(/3)|(3a 2 -4a) = -V59 » 3.84 

of the fundamental domain of P. And because of the periodic tiling of the plane with Q as 
a prototile, it is also easy to calculate the area // int (fi) of f2: as /i int (9f2) = 0, /i int (f2) equals 
the area of a fundamental domain of the corresponding lattice of periods. This gives 

(22) Mint (fi) = |det((P 2 -P 6 ),(P 2 -P 3 ))| = |Im(/3)|(a 2 -a) = iJL(3a 2 -2a+17) ra L77 " 

^ V59 

Then the following lemma applies. 

Lemma 1. Let P be a lattice in M x IR m , |P| 6e the volume of a measurable fundamental 
domain of P in I x M m u>ii/i respect to the product [i ® fx int o/ the Lebesgue measures jx, 
fi int on M, M m , respectively. Assume that we have a cut-and-project scheme like in Q. If 
Q is a bounded subset of M m with almost no boundary, then the density dens(A(0)) of the 



7 Note that the discriminant of Q(a) is —59. The volume \r\ is proportional to the square root of the 
absolute value of the discriminant. The proportional constant is one factor of | because there is one complex 
conjugate pair (3,/3 of algebraic conjugates of a, see |1 II Chapter II, Section 4.2, Theorem 2]. Note that we 
also have a formula for | Im(/3)| in terms of a by this: 

|Im(/3)| = 3_(-8a 2 + 25a-6). 

2-^59 
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corresponding regular model set in R is 

^int(^) 



dens(yl(n)) 



Proof. This follows from [28, Proposition 2.1] because the projection ir is one-to-one on r by 
construction. □ 

With (|H|l. (|21|) and (|22jl. it is now easy to check that the density of the model set A(Q) and 
the density of SKol(3, 1) are equal, i.e., 

f0*\ ^int(^) _ 1 

Proposition 6. The sequence SKol(3, 1) is a subset of A{Q). Further, they differ at most 
on positions of zero density and therefore have the same pure point diffraction spectrum. 

Proof. We choose 0, —v B € r. Then their projections into internal space are elements of the 
attractor O, because ^(0) = and / 3 (— /3) = — 0. But starting with these two points, the 
iteration of the IFS in internal space just corresponds to the iteration in (J3J), respectively 
in physical space. Therefore, SKol(3, 1) C A(Q), because the star map of all iterates of 
and —a (i.e., and — /3) stay in 0. Equation ([2*3*)) shows that both sequences have the same 
density. So they can at most differ on positions of zero density. 

Regular model sets have a pure point diffraction spectrum, see [*"] El 1239 an d references 
therein. Therefore, the diffraction spectrum of A(Q) is pure point, and SKol(3, 1), differing 
at most at positions of zero density, has the same spectrum by an argument in |18j . □ 

Theorem 1. EKol(3, 1) is a regular model set (except possibly for positions of zero density) 
and has a pure point diffraction spectrum. Its autocorrelation is a norm almost periodic point 
measure, supported on a uniformly discrete subset ofX[a\. 

Proof. The first assertion follows from Proposition El The autocorrelation measure (see |18| 
for details) is supported on 71(17) — A(£l) C Z[a]. Since A(U) is a model set, it is also a Meyer 
set, hence A(Jl) — A(O) is Delone. The norm almost periodicity follows from 6, Theorem 
5]. □ 

What we have proved so far is enough to calculate the diffraction spectrum of SKol(3, 1) 
and Kol(3, 1), see Section"""! But in the next section, we want to show that A(Q) really equals 
£Kol(3,l). 

Remarks: For unimodular substitutions of Pisot type, i.e., Kol(2m± 1, 2m =p 1) with m > 1, 
the procedure is essentially the same. Unfortunately, the IFS does not decouple like in (jl2j) for 
m > 1, which makes it technically more involved. For non-unimodular substitutions of Pisot 
type, the internal space is more complicated in having additional p-adic type components, 
see [23 El HI US] for further details and examples. 

The sequences £Kol(p, q), which are not of Pisot type, do not have a pure point spectral 
component outside k = by an argument in pil 5 see |31j for details. 

Some of the results given have been studied extensively under the name of "Rauzy fractal", 
e.g., that the windows have non-empty interior, that the windows do not overlap in this case 
(this follows from the so-called strong coincidence condition) and also the periodic tilability 
seems to follow from results in [TT 1121 13*5] . But we also need the lattice of the periodic tiling 
explicitly, as well as the induced IFS (|16|) for the boundary. Therefore, we opted to give an 
elementary and complete derivation here. 
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4. £Kol(3, 1) IS A GENERIC MODEL SET 

o o o o o o 

We set f2 = £l A U £1 B U Q, c C £1 and 90 = O \ O = dQ A U <9r2 B U d£l c . Then we can improve 
a statement of Proposition El 

o 

Proposition 7. £Kol(3, 1) is equal to the model set A(£l). 

o 

Proof. £Kol(3, 1) C A(Q): Note that the mappings f i (i G {0,1,3}) of (JTUt ) are similarities 
(all directions are contracted by the same factor, here \(3\). Therefore, they map balls around 

o o 

x to balls around f^x). Furthermore, they map balls in £l k to balls in Cl i (k,£ G {A, B,C}). 
Since the starting points of the iteration Q in internal space, namely —(3 and 0, are inner 
points of Q B and fl A by Corollary one can also find balls of radius e > around —j3 

o o 

and which lie entirely in Q B and £l A , respectively. Since the iteration in physical space 
corresponds to the IFS in internal space, the star map of an arbitrary point in SKol(3, 1) is 

o 

thus a point of fi. 

A(h) C SKol(3, 1): Suppose x G yl(0)\EKol(3, 1). Then x* and all its iterates of the IFS (TTU1) 

o 

are in f2, by the same reasoning as before. Furthermore, the mappings f i (i G {0,1,3}) are 
affine similarities and therefore all iterates of x* are disjoint to all of the iterates of —(5 and 
0. But then 

densSKol(3, 1) < densyl(O), 

o 

because the set of iterates of x under inflation has positive density in A(Vt). This contradicts 
Proposition El □ 

Note that not only the original sequence with 3's and l's is inversion symmetric (see @), 
but also the positions in A(Q AB ). 

Corollary 3. A(U AB ) is inversion symmetric in the center — \ct. 

Proof. The starting points —(3 and and the IFS H12|) are inversion symmetric in the center 
P 5 . This corresponds, in physical space (by Galois conjugation), to inversion symmetry in 
the center —\a. □ 

For the following, we need some more definitions, see (Hj and |26| . If A is a discrete point 
set in M. d , we call V r {u) an r -patch of a point u G A, if V r (u) = A D B r (u), where B r {u) is 
the ball of radius r about u. Often, we are only interested in the set of {V r {u) \ r > 0, u G A} 
and an element of this set is simply called a patch. Two structures A 1 and A 2 are locally 
indistinguishable (or locally isomorphic or LI) if each patch of A 1 is, up to translation, also 
a patch of A 2 and vice versa. The corresponding equivalence class is called Ll-class. 

A discrete structure A is repetitive, if for every r > there is a radius R(r) > such 
that within each ball of radius R(r), no matter its position in there is at least one trans- 
late of each r-patch. Note that every primitive substitution generates a repetitive sequence, 
wherefore SKol(3, 1) is repetitive. 

We now look at generic model sets and show that A(Q) is actually generic. 

o 

Lemma 2. The set C = {cGK 2 |(c + dQ) n L* = 0} is dense in R 2 , especially c + L* C C 
for c G C. 
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Proof. The set C is not the empty set by Baire's category theorem (d£l is a meager set, L* is 
countable), by standard arguments, which in this context first appeared in |27| Section 2.2.2], 
also see But if c G C, then c + 1 G C, Vt £ L* (notice that L* is an Abelian group), and 
L* is dense. □ 

Proposition 8. 

(i) The model set A(c + fi) is repetitive and generic for c £ C. 

(ii) T/ie model sets A(c + O) and A{c + O) are LI for c,c G C. 

(iii) EKol(3, 1) and A(c + fij are LI for c£C. 



Proof, (i) The model set A(c + fi) is generic by the definition of the set C. It is repetitive 



by |28[ Theorem 6] and |29[ Proposition 3.1]. 



(ii) This is a by now standard argument, apparently first used in |27[ Lemma 2.1]. 

(iii) Since A(c + 0,) is repetitive, by |27| Lemma 1.2] it is enough to check that every patch of 

o 

SKol(3, 1) also occurs as a patch of A(c + fi): Let V be a patch of EKol(3, 1). Then V* C £1 
by Proposition [7[ and since V is a finite patch, we even know that there is an e > such 

that V* C O e , where fi £ = {i G M 2 | t G CI, dist(t,dft) > e}. The set C is dense by LemmaEJ 
therefore there is a c G C such that V* C c + fl. Then P C A(c +J2) and EKol(3, 1) and 



vl(c + J7) are LI. Since LI is an equivalence relation, it follows from (ii) that SKol(3, 1) and 
A(c + Q) are LI for every c G C. □ 

Proposition 9. Define C = {cSC|OGc + O}. Then, for every r > and c G C , we get 

b p (o) n SKoi(3, i) = £ r (o) n (yi(c + o) + 1) = b p (o) n A(c + t* + n) 

for an appropriate t G Z[a]. Additionally, it follows that c+t* G C . 

Proof. By Proposition [ ^ iii)] we know that A(c+Cl) and SKol(3, 1) are LI, therefore the patch 
B r (0) n EKol(3, 1) occurs somewhere in A(c + fi). By the choice of C , we can translate this 
patch in yl(c + Q) with i G Z[a] to the origin. In internal space, this is a translation t* and 
by Lemma |21 we have c + 1* G C. But EKol(3, 1) has a point at the origin, so we even have 
c + t*£C . □ 

Theorem 2. EKol(3, 1) is a regular generic model set. 

Proof. Since EKol(3, 1) is the fixed point of a primitive substitution, it is repetitive, and the 
corresponding dynamical system is minimal, see |29[ Proposition 3.1]. 

By Proposition P ^iii)! the model set A(c + Q) for c G C is in the Ll-class of EKol(3, 1), 
hence the latter must be the limit of some sequence (t i ) i of translations of A{c + Q), where 
the translates t i can be restricted to elements of Z[a] and c + t i G C by Proposition El 

O o 

However, if <9f2 fl L* ^ 0, there is a point x G dQ nL*nfi with x G dQ A n dQ B (so x lies 
on the common boundary of £1 A and Q B ). This is because by appropriate combinations of 

the mappings g 1 ,g 2 ,g 3 ,r, K, f , f v f 2 , f 3 , f 4 and the translations by P 2 - P 1 , P 2 - P 3 , P 2 - P 6 

o 

of Sectional which all map L* onto L*, we can "move" points on dfl from every "edge" to 

o 

every other "edge" 9 . We have e = dist(x, dU \ (dQ A n dQ, B )) > 0. 

^That every patch of A(c + tl) is also one of EKol(3, 1), follows then together with the repetivity of A(c + Q,). 
q o ' o 

We even get that, with one point x £ dQ n L , there is a dense set of points in dfl n I/*, because we can 

always "move" x to the edge [P2>-P3]i a PPly the IFS 1161 there and "move" this edge, with now dense points, 

to every other edge. 
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The inverse star image of this point must then be in any limit of sequences A(t* + c + fi) 
with c + t* —* 0, but it is not in EKol(3, 1) — which is a contradiction. So no such point can 
exist and d£l fl L* = 0. 

Regularity was established in Theorem ^ together with Proposition [Tj □ 

This argument is rather general and applies in other situations as well. For a more elemen- 
tary proof see the appendix. 

Remark: By Proposition we know that (JteG^ + ^) = ^ 2 > where G = (—0 + 1, /3 2 — 2(3}% 
is a rank 2 free Abelian group (a 2-dimensional lattice) and by Theorem |2] that L* = Z\fl] C 

o 

Utec(* + w here U denotes disjoint union. In physical space, we get 

L = Z[a] = {J (s + SKol(3, 1)), 

where 

G' = (-a + 1, a 2 - 2 a) z = (£ c - £ B , t A - £ B ) Z , 

i.e., Z[a] is the disjoint union of translates of the regular, generic model set EKol(3, 1). The 
set of translations needed is a rank 2 subgroup, whose Galois dual in Z[/3] is a lattice. But 
we can also write 

L = Z[a}= (j (r + C), 

reEKol(3,l) 

so L/G' is a coset system with the structure of a model set. Now, let A(m) be the m-th 
element of SKol(3, 1) (m G Z). Then one can show that the induced group structure on this 
coset system is (A(m)) + (A(n)) = (A(m + n)), i.e., it is the action of Z on EKol(3, 1). This 
group structure lines up with the deformation in the next section. 

5. Deformation and Diffraction 

In the cut-and-project scheme (J7J) with H = M m , let (p : W 11 — ► M. d be a continuous function 
with compact support (e.g., O). We call 

A v (n) = {x + <p(x*)\x€A(n)} 

a deformed model set if it is also a Delone set, see The model set A(£l) can be seen as 
deformed model set where the associated function ip is trivial, i.e., (p = 0. The diffraction 
spectrum of (deformed) model sets (where each of its points is represented by a normalized 
Dirac measure, say) can be calculated explicitly, see jH] for details. We write 5 k for the Dirac 
measure at k, i.e., S k (f) = f(k) for / continuous. Also, we need the dual of a lattice r C R n 
defined as 

r* := {y G R n | z.y G Z, Vx G T}, 
with x.y denoting the Euclidean scalar product. 

Proposition 10. jH] Let A (Q) be a deformed model set in M rf constructed with a regular 
model set A(Q) and a continuous function ip of compact support. Then, the diffraction pattern 
of A (Q) is the positive pure point measure 

7= \c k (A^(n))\ 2 8 k , 

ken(r*) 
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where r* is the dual lattice of T, 5 k is the Dirac measure at k and c k (A (Q)) is the Fourier- 
Bohr coefficient of A^Cl) at k. This Fourier-Bohr coefficient exists and has the value 

f If p -2wi(k.ip(y)-k*.y) J. if ( k F) f f* 

(24) c k (A v (Q)) = \ n ^ ^ V^ } k)€l , 

I U, otherwise. 

□ 

For a regular model set A(£l) (where tp = 0), the Fourier-Bohr coefficient is just given by the 
(inverse) Fourier transform of the characteristic function of the window Q. 

For £Kol(3, 1), we note that the support F of the spectrum is dense in R since it is given 
by the Z-span of the projection of the dual lattice vectors, i.e., 

(25) F = Z-tt(v* a )+Z-tt(v* b )+Z-'k(v* c ), 

but it{v* B ) = {a — 1) n(v A ) and therefore they are linearly independent over Q. 

To deform SKol(3, 1) to Kol(3, 1), we make the linear ansatz <p(x*) = ax\ + bx^, where 
x* denotes the ith Cartesian component of the vector x* G M 2 . With this <p, we now deform 
all bond lengths t i to the average bond length £, i.e., we have to solve the following linear 
system of equations (i £ {A, B, C}): 

+ a (vi) 2 + b (vi) 3 = £. 
This over-determined system is solved by 

a = £ - 1 = 2 (-a 2 + a + 5) rs 1.17 and 



Im([3) , , ,1 

— — -(-a 2 -17a + 31 = 410 ~ 2 

59 * 



(-a 2 - 17a + 31) = T (-413 + 885 a - 59) » -0.13. 

2^59 



Due to the linearity of y> (and the positivity of the bond lengths involved), this deformation 
does not alter the order of the points (i.e., for x, x' £ A(£l) with x < x', we always have 
x + tp(x*) < x' + ip{x' )). Note that the support F of the spectrum stays the same as in 1)25)1 . 
only the Fourier-Bohr coefficients change. 

The positions in vl (QJ are now subsets of £ • Z. To be more precise, we even have 

A^n A )uA v (n B )uA v (n c ) =£■%. 

Because of this embedding into £ ■ Z, the diffraction spectrum of each of the aperiodic sets 
ApiClj) is (£ ■ Z)*-periodic [2], i.e., it is periodic with period l/£ (note that Z* = Z; the 
diffraction spectrum of SKol(3, 1) is not periodic). This might not be obvious from 1)24)1 at 
first sight, but for n £ Z we have (note that 7r(i>*) = pjt) 
n 

(26) - = nir(v* A + v* B + v* c ) = ir(nv* A + nv* B + nv* c ), 

therefore, for every k' = k + j with (k, k*) £ r*, there is also a k' with {k' , k' ) £ r* given 
by 

k'* = k* + 7r int (n v* A + nv* B +nv* c ). 
But with the chosen <p we get 



(27) (k' - k).<p(y) - {hT - k*).y 



ny x 



a 



{■K\nt{v*A + V B + V* c )) 1 



ny 2 



J ~ (7Tint(«A + V *B +V* c ))r 
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because each of the terms in square brackets vanishes. Therefore 

holds, and the spectrum is periodic with period \ 

To obtain the diffraction spectrum of Kol(3, 1) from here, one only has to rescale the 
positions in A (£1) by a factor of l/£. To summarize: 

Theorem 3. The bi-infinite sequence Kol(3, 1), represented with equal bond lengths, is a 
deformed model set and has a pure point diffraction spectrum. □ 

Remarks: By the same method, we can also find a deformation <p(x*) = ax* + bx\ such 
that we represent the letter T' of Kol(3, 1) with an interval of length i and the letter '3' with 
one of length 3£. For this, the letters A, B, C have bond lengths 61, At, 2t, respectively. For 
the parameters of the deformation (the average bond length must be I again), we get 



Tm(J3) 
59 ' 



15 a + 11 



0.49, 



15 a- 1] 



-0.016 and 



(-179 a 2 + 379 a + 3) 



2V59 



: (1239 a 2 - 767 a -4661) 



-0.36. 



Now ((2*6^1 changes to 



71 



— = ir(6nv* A + Anv* B + 2nv* c ) 



and with the same calculation as before one gets an equation which corresponds to Q27j) , where 
the two terms in square brackets 



j-(^v* A + 4v* B + 2v* c )) 1 



and 



~-(TT int (6v A + Av B + 2v* c )) 2 



also both vanish. Therefore, the spectrum is periodic with period l/l as expected j2], since 
Afc(Q) ^ £ ■ 7L. This representation with integer bond lengths (after rescaling) has the advan- 
tage that the union of the three aperiodic sets A^i^l^) is still an aperiodic set. Clearly, it is 
also pure point diffractive. 

Kol(3, 1) in its natural setting with intervals of length 1, or of lengths 3 and 1, can be 
obtained as a deformation of the model set SKol(3, 1) derived above, where the intervals 
have incommensurate length. The basic theory of this is fully developed in ||JE3, but one 
can also understand, from a dynamical systems point of view, which deformations are stable 
in the sense that they do not change the spectral type of the dynamical spectrum (and hence 
of the diffraction spectrum, due to unique ergodicity), see jl]. 
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Appendix: An alternative proof of Theorem [2 

By Proposition El we can choose a sequence (c i ) i with q = c + t* G C (i* G L*) such that 
£? r n EKol(3, 1) = B r n A[c i + 0) for every sequence (r i ) i with > i. Also, this statement 
holds for every subsequence (q . ) j . 

O O o 

Now, assume dQnL* ^ 0. Then we have a point x G <9r2nL*nSl with x G cT^ndO^. Set 
e = dist(x, dQ \ {dn A n 5fJ B )) > 0. Then a translation y / of Q with y G B e (0) n C has 

o 

the following effect: x G y + fi, because by the definition of C, x cannot be on the boundary 

o 

y + d£l, and by the choice of e , it must either be in y + £l A or y + 

Now take a sequence (cjj as above. Clearly, this sequence must converge to 0. Therefore, 
there is an N such that |cj < e for all i > N. By choosing an appropriate subsequence (c^.)^ 
we get a sequence (c^ with = c { . such that x is always either in &j + A or in c- + B . 
Also we have B r nSKol(3, 1) = B r nA(c i + Q) for r i > i. But both ir and 7r int are one-to-one. 
Therefore, the inverse of the star map of x must be a point of each A{c i + 0) and it also must 
be in B R (0) for some R < oo. But by Proposition [7J it is not in SKol(3, 1). Therefore we 
get a contradiction and our assumption is wrong. So, G C and SKol(3, 1) is generic by 
Proposition | ^i)| □ 
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